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Mesoscopic thermoelectric heat engine is much anticipated as a device that allows ns to utilize 
with high efficiency wasted heat inaccessible by conventional heat engines. However, the derivation 
of the heat current in this engine seems to be either not general or described too briefly, even 
inappropriate in some cases. In this paper, we give a clear-cut derivation of the heat cnrrent of 
the engine with suitable assumptions beyond the linear-response regime. It resolves the confusion 
in the definition of the heat current in the linear-response regime. After verifying that we can 
construct the same formalism as that of the cyclic engine, we find the following two interesting 
results within the Landauer-Biittiker formalism: the efficiency of the mesoscopic thermoelectric 
engine reaches the Carnot efficiency if and only if the transmission probability is finite at a specific 
energy and zero otherwise; the unitarity of the transmission probability guarantees the second law 
of thermodynamics, invalidating Benenti et al.’s argument in the linear-response regime that one 
could obtain a finite power with the Carnot efficiency under a broken time-reversal symmetry. These 
results demonstrate how quantum mechanics constrains thermodynamics. 

PACS numbers: 05.60.Gg, 73.23.-b, 05.70.Ln, 72.15.Jf 


I. INTRODUCTION 

Thermoelectric heat engine [M] is much anticipated 
as a device that allows us to utilize wasted heat inacces¬ 
sible by conventional heat engines. This engine operates 
in a nonequilibrium steady state and converts heat to 
useful electrical power steadily, so that we do not need 
nonsteady processes used in cyclic engines, such as adi¬ 
abatic compression, isothermal expansion, and so forth. 
Its efficiency, however, so far has been too low to use 
in terms of the figure of merit Z{T), which is a serious 
problem in the field of thermoelectricity [5TH55] . 

The mesoscopic thermoelectric heat engine [DI m- 
ini IB ITOHTOl m]-^ has emerged as a possible solution. 
This engine is expected to have high efficiency thanks 
to the potential of nanoscale thermoelectricity [?TM33j . 
Moreover, this engine can be a powerful tool to inves¬ 
tigate how quantum mechanics affects thermodynamics. 
For example, it has been argued that the heat current 
may be bounded because of the uncertainty principle 
[IB12111B1M] and that the unitarity of the scattering 
matrix may give a new bound for the Onsager coefficients 
in the linear-response regime [niiiaES]. 

In order to calculate the efficiency of a mesoscopic 
heat engine, we need to know the expression of the 
heat current going into or from a mesoscopic system, 
such as a quantum dot and a quantum wire. Although 
many researchers have already used the definition of 
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the heat current in the linear-response regime |1Q| I13h 
[iiiiiiiiaisBED] as well as in nonlinear-response regimes 

[31 SHi IB im El EMB EB El EBEB ESHU, they 

did not give its derivation or explanation in most of the 
papers. There are several papers in which they explain 
the heat current. In Ref. |42j . for example, Sivan and 
Imry gave a kind of derivation, which is the only one in 
the linear-response regime as far as we know. However, 
it is not clear from the present authors’ point of view; 
the derived expression is also inappropriate in nonlinear- 
response regimes as we will show below. In nonlinear- 
response regimes, although several authors described the 
heat current [28l EB 1381141] , their explanations seem too 
brief to understand for readers who are not familiar with 
the mesoscopic heat engine. Moreover, since the presen¬ 
tation was done mostly in the Landauer-Biittiker formal¬ 
ism, the general framework seems to be hidden. 

We, in the present paper, construct a thermodynamic 
formalism beyond the linear-response regime under rea¬ 
sonable assumptions. Since we do not use the Landauer- 
Biittiker formalism in the first half, our general formu¬ 
lation is applicable to systems with interactions as long 
as the assumptions are satisfied. In Sec. II, we explain 
that thermoelectricity works as a heat engine under suit¬ 
able conditions. In Sec. IH, after we briefly overview 
the confusion in the definition of the heat current in the 
linear-response regime, we first show the derivation of the 
heat currents in the thermoelectric steady-state heat en¬ 
gine with suitable assumptions. With the heat currents 
that we derived, we construct the general formalism. In 
Sec. IV, we give an example of the heat engine using the 
Landauer-Biittiker formalism, which we call the meso¬ 
scopic thermoelectric heat engine. We can confirm the 
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FIG. 1. Schematic picture of a mesoscopic heat engine. We 
set the chemical potential of the right reservoir higher than 
the left, while the temperature of the left reservoir higher than 
the right so that an electric current may go from left to right 
against the difference of the chemical potential. 


nonnegativity of the entropy production of this engine. 
Moreover, we find that the efficiency of the mesoscopic 
thermoelectric engine reaches the Carnot efficiency if and 
only if the transmission probability is finite at a certain 
energy and zero otherwise. In Sec. V, we give a model 
that expresses the situation when inelastic scatterings oc¬ 
cur in the central system in Fig.[l] In Sec. VI, we consider 
the system with a broken time-reversal symmetry using 
the Landauer-Biittiker formalism. We find in nonlinear- 
response regimes that the unitarity of the transmission 
probability invalidates Benenti et al.’s argument [5] in 
the linear-response regime that one could obtain a finite 
power with the Carnot efficiency under a broken time- 
reversal symmetry. It is remarkable that the unitarity, a 
quantum-mechanical concept, constrains thermodynam¬ 
ics; it guarantees that the Carnot efficiency is achieved 
only at zero power. 


II. THERMOELECTRIC HEAT ENGINE 

Let us explain the thermoelectric steady-state heat en¬ 
gine. Consider a central system, for example quantum 
dots or quantum wires, attached to two reservoirs on 
both sides; see Fig. We then make the following three 
assumptions: (i) the reservoirs are so much larger than 
the central system that it is always at equilibrium even 
if they interact with the central system, and hence we 
can define thermodynamic quantities of each reservoir, 
such as the temperature and the chemical potential; (ii) 
the central system has reached a nonequilibrium steady 
state in which there are constant flows; (iii) there is no 
entropy production in the central system because elec¬ 
trons undergo only elastic scatterings there. 

We can regard this system as a heat engine in the fol¬ 
lowing situation. We set the chemical potential of the 
right reservoir higher than that of the left, while the 
temperature of the left reservoir higher than that of the 
right. We particularly set the reservoirs as well as the 
central system so that an electric current may go from 


left to right against the difference of the chemical po¬ 
tential. What happens per unit time is the following. 
Electrons gain heat Jq from the hot left reservoir, flow 
to the right against the potential difference, during which 
electrons do the work of amount 

W = IV, (1) 


where 


II 

(2) 

is the electric current and 


MR - Ml 

( 3 ) 


is the voltage difference, and then dump heat Jq to the 
cold right reservoir. We can thus consider this system as 
a heat engine. This is an explanation specific to the case 
of electrons, which we can make more general as follows. 
The central system gains heat from the left reservoir, does 
work, then dumps heat to the right reservoir. In this per¬ 
spective, we can regard the central system as a working 
system of cyclic heat engines. We will use this perspec¬ 
tive hereafter. Because there is no entropy production in 
the central system, we must have 

4 --^^ = ( 4 ) 

Its efficiency tj is thereby given by 


W 

JL =1- 


4 

4 ’ 


( 5 ) 


which is the same as the standard cyclic heat engine. 

Note that we can regard the system as a heat engine 
only when Jjv and Jq are positive. The electrons do not 
necessarily flow from the hot reservoir to the cold one; 
the current which goes to the left is considered to be 
negative. The direction of the flow depends on ml, Mr, 
Tl, and Tr as well as details of the system. For example, 
we can regard the system as a refrigerator when electrons 
go from right to left against the temperature difference. 
In this case, the efficiency, which is called the coefficient 
of performance, is given by Mcop = 4/^^^^^ ™ 

Eq. 0. 


HI. HEAT CURRENT OF THE 
THERMOELECTRIC HEAT ENGINE 


A. Confusion in the definition of the heat current 
in the linear-response regime 

Based on the above argument, we realize that we need 
two heat currents Jq and Jq in order to discuss the ther¬ 
moelectric heat engine properly. Before deriving these 
heat currents, we mention the confusion in the definition 
of the heat current in the linear-response regime. 
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We first note that the energy current Je is often re¬ 
ferred to as a “heat” current [431 Si] ■ This would be cor¬ 
rect, though confusing, if the electrons did not do work 
and hence all energy became heat. This is certainly not 
correct in the situation in Fig. 

Another definition Jq = Je — ^J■JN was often used in 
the dawn of the research of the heat current in mesoscopic 
systems [12] • This definition may have been taken from 
Eq. (17.8) in Callen’s textbook [3^. Since this “heat” 
current was not microscopically derived, we do not clearly 
know where it flows. It is indeed ambiguous of which part 
of the system in Fig. I^the chemical potential /r of the heat 
current Jq = Je — ^^JN is. We should probably choose /i 
so that Jq may satisfy Onsager’s reciprocal theorem. For 
example, in Ref. [42] the authors chose /x as (^l + Mr)/ 2 
and in Ref. [33] the author chose /x as /xl. The choices 
do not make any difference in the linear response of the 
voltage difference V = (^r — between Jq and Jq 

but differ in higher orders. 


B. Derivation of the heat currents going into or 
from a reservoir 

We here derive the heat currents Jq and Jq thermody¬ 
namically consistently. Although many researchers have 
used them, the derivation is often skipped or described 
too briefly. 

Our starting point is to identify the heat currents com¬ 
ing into or going out of the central system as those coming 
from or going into the reservoirs. We will give an impor¬ 
tant remark on this point later in this section. We now 
derive the latter using general thermodynamics. The first 
law of thermodynamics gives the relation 

dC/“ = dQ“-bdW“, (6) 

where a = L,R, dt/“ and d(5“ denote the energy and 
heat flowing into the left or right reservoir, and dIF“ = 
^o,dA^“ is the work done to the reservoir. Using Eq. Q, 
we can express d(5“ in the form 

dg“ = dU“-AiadAf“. (7) 

Since we treat the nonequilibrium steady state, we can 
define the changes of the particle number and the energy 
in a reservoir as a steady current, which enables us to 
define the particle and energy currents as 

^ dN^ dN^ ^ dU^ dU^ 

where we used the conservation of particle number and 
energy, dA^^ + dA^^ = 0 and dU^ + dU^ — 0, because 
there is no dissipation in the central system. The nega¬ 
tive signs appeared because we defined the positive direc¬ 
tion so that currents going to the right may be positive. 
We then define the heat currents using Eqs. Q and ^ 


in the form 

Jq = -^ = Je - fiL Jjv, 

dO^ 

Jq = = Je - ^Jn, (9) 

where Jq is the heat current flowing from the left reser¬ 
voir into the system and Jq is that flowing from the 
system into the right reservoir. 

Note that these expressions of the heat currents are 
valid in the presence of many-body interactions if our 
assumptions (i) to (iii) hold. We can therefore apply 
them to systems for which we cannot use the Landauer- 
Biittiker formula. We will also discuss a possible formu¬ 
lation in Sec. V when there is an entropy production in 
the central system. 

Let us come back to our starting point above. It is 
crucial to note that the working system of the engine is 
now the central system in Fig.j^ for which we cannot de¬ 
fine the temperature and the chemical potential because 
it is highly nonequilibrium. In the textbook [37], Datta 
argued that we could define an “effective” chemical po¬ 
tential in quantum wires, but it is in fact not a thermody¬ 
namic observable. We should therefore be aware that it 
is not trivial at all to define the heat and the work for the 
central system. We here identify the heat current coming 
into the central system as that coming from a reservoir 
under the assumption below, which corresponds to the 
one that there is no entropy production in the central 
system. 

According to the first law of thermodynamics, in or¬ 
der to calculate the heat current coming into or from 
the central system, we need to specify the energy current 
coming into or from it and the work done to or by it. 
The former is easy to specify as Je because the energy 
current is a conserved quantity. The problem arises when 
we try to specify the work done to or by the central sys¬ 
tem, in which we cannot define thermodynamic intensive 
variables. 

We here make the following assumption to specify it: 
the work done to or by the central system is equal to 
that done by or to the reservoir. This assumption lets 
us find that the work done to the system on the left side 
is equal to that done by the left reservoir, /xlJjv, and 
the work done by the system on the right side is equal 
to that done to the right reservoir, /xrJjv. We can thus 
find that the expression of the heat current coming into 
or from the central system is equal to those that we de¬ 
rived in Eq. Although the expression of the heat 

currents corresponds to that used in previous researches, 
we believe that our derivation is more accurate than the 
previous ones. The crucial point is that we derived the 
heat current coming into or from the central system, not 
a reservoir; we need to make the assumption above to 
specify it. The assumption is probably equivalent to the 
assumption (iii) in Sec. 11 that there is no entropy pro¬ 
duction in the central system. 
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C. General formalism of the thermoelectric heat 
engine 


Let us show that we can construct the same formal¬ 
ism as that of the cyclic engine using the heat currents 
Eq. We first easily find that they indeed satisfy 

Eq. Q. This is the first check of the consistency of the 
expression Eq. (|^. 

We then show that the upper limit of the efficiency 
Eq. (§ is the Carnot efficiency as is expected from the 
theory of the standard cyclic heat engine. Let 

^Ldt = ^ (10) 

and 

Sndt = ^ ( 11 ) 


denote the entropy productions in the left and right reser¬ 
voirs, respectively. Using Eq. § , we can relate these 
entropy productions to the heat currents as 

Jq = -Tl5l, 4 = TnSn- (12) 

According to the second law of thermodynamics, an en¬ 
tropy production of an isolated system increases. Our 
whole system, which consists of the two reservoirs and 
the central system is indeed isolated. The net entropy 
production of the whole system 

S = Si^ + Sk (13) 


is thus non-negative, that is, —S'r < S'l; we assumed that 
there is no entropy production in the central system. We 
thus have from Eqs. Q, ([^, and (12), 


V = 


- TnSn 


—T’l*5l 


^ 1 


(14) 


where rjc is the Carnot efficiency. We can achieve the 
equality if and only if = — «5r, that is, 5 = 0. We 
conclude that with the heat currents Eq. § , we can pro¬ 
duce Eqs. 0 and ( [T4| , which are the same formulas as 
those of the cyclic engine. 


IV. HEAT CURRENT IN THE MESOSCOPIC 
HEAT ENGINE USING THE 
LANDAUER-BUTTIKER FORMULA 

In order to justify the definitions Eq. ^ further, we 
here derive microscopic expressions for the mesoscopic 
thermoelectric heat engine. We additionally assume here 
that the central system in Fig.[^is a noninteracting coher¬ 
ent conductor that accommodates the Landauer-Biittiker 
formalism m- Note again that the arguments in Secs, 
n and HI are valid even in the presence of many-body 


interactions, to which the Landauer-Biittiker formalism 
is not applicable. 

We can obtain and Je in Eq. 0 microscopically 
as 

^ y *'r(e)(/L(e) -/R(e)), (15) 

JE = ^j de T(e)e(/L(e) - /R(e)), (16) 

where h is the Planck constant, 


I -I- exp(/3a(e - Ha)) 


(17) 


is the Fermi distribution function of the reservoir (a = 
L,R), Pa = T~^ is the inverse temperature, and r(e) is 
the transmission probability at energy e. Substituting 
Eqs. (15) and (16) into Eq. (|^, we arrive at 


’^Q = \ J “ Ma)(/L(e) - /R(e)), 


(18) 


where a = L, R. 

Note that these heat currents satisfy Onsager’s recip¬ 
rocal theorem when we expand them in terms of ap¬ 
propriate affinities; we can verify the theorem by ex¬ 
panding Jn and Jq in terms of = /3 l(ml — Mr) 
and = -(/ 3 l - Pk) or Jn and in terms of 

= ^r(ml - Mr) and = -(/ 3 l - ^r) gS]. 

We can also verify the non-negativity of the entropy 
production with the heat currents Eq. 0 and their mi¬ 
croscopic expressions Eq. (18), although it must be sat¬ 


isfied anyway according to the second law of thermody¬ 
namics. Using Eqs. (12) and (13), the total entropy pro¬ 


duction of the system is generally expressed with the heat 
currents: 


• • • 4 Jq 

S = Si. + Sn = -T^ + ^. 

J-L Jr 


(19) 


Substituting the microscopic expression of the heat cur¬ 
rents Eq. (18), we obtain the microscopic expression of 


the total entropy production in the form 


5 = 


de T(e)(/L(e) - /R(e))log 


/L(e)(l - /R(e)) 


./R(e)(l-/L(e))J’ 

(20) 

which Whitney derived |18j starting from the nonlinear 
Landauer-Biittiker formalism. Since r(e) > 0, the in¬ 
tegrand in Eq. (20) is always nonnegative, which leads 
to the nonnegativity of the entropy production. The 
fact that we reproduced Whitney’s microscopic expres¬ 
sion Eq. ( |20| ) endorses our general expressions of the heat 
currents Eq. 0. 

Since we know that we can achieve the Carnot effi¬ 
ciency when the total entropy production of the whole 
system is zero, let us find the transmission probability 
that satisfies 5 = 0. We easily see the following condi¬ 
tion: for each value of e, T(e) = 0 or fh{e) — fnie) = 0. If 
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FIG. 2. The energy filter transmission probability. 


FIG. 3. A three-terminal model. The third bosonic reservoir 
on the top represents the energy dissipation. 


r(e) = 0 for any e or if fh{i) — /R(e) = 0 for any e, how¬ 
ever, the transport would not happen, the engine would 
not work, and the efficiency rj would be trivially 0/0. 
The only nontrivial condition is given by T(e) 0 and 
— /r/e) = 0 at one value of energy; if we demanded 
/L(e) = /u/e) at two values of energy, they would be 
equal at any values of energy. 

In order to achieve the above condition, the transmis¬ 
sion probability should be of the form 


we here set c in Eq. (211 as unity for simplicity. The 
power is given by [55] 


— Ml) = 


Mr ~ Ml 
h 

Mr ~ Ml 
2h 


r<^c+A 


/ de{fi^{e) - fnie)) 

ec 

E(ec)A2 + 0(A3), (24) 


r(e) 


c for e = Cc 
0 otherwise, 


( 21 ) 


where c is a constant satisfying 0 < c < 1 and 
TlMR ~ TrMl 


Cc = 


Tt. - Tr 


( 22 ) 


is derived from the condition fhi^c) = /R(ec)- Using 
Eqs. ([^, (|^, (15), and (18), we indeed find that this 
transmission probability is a necessary and sufhcient con¬ 
dition for the Carnot efficiency: 


Ml - Mr 
ec - ML 



= Me- 


(23) 


The transmission probability Eq. mi) for the Carnot 
efficiency that we derived is well known in the litera¬ 
ture on thermoelectricity [niiiiiniiiTiiiiiMiiin]- How¬ 
ever, in the mesoscopic thermoelectric engine with the 
Landauer-Biittiker formalism that we treat in the present 
paper, there has been no discussion that the transmis¬ 
sion probability Eq. ( |21[ ) is a sufficient condition for the 
Carnot efficiency [UlTa [Ml HH], as far as we know. In 
contrast, we showed that the transmission probability 
Eq. (21) is a necessary and sufficient condition for the 
Carnot efficiency straightforwardly from the expression 
of the entropy production Eq. (20) along with the condi¬ 
tion 5 = 0. 


In order to see what limit of this engine corresponds 
to the quasistatic limit of the cyclic heat engine, let us 
consider the transmission probability Eq. mi in the limit 
of A —>■ 0 of the energy window [Ml HH] shown in Fig. 


where F{e) = 9/i9e(/L(e) — /R(e)). The heat current 
Jq = (cc — ML).C(ec)A^/2M -I- 0{A^) vanishes and the 
efficiency rj = rje — C’(A) becomes the Carnot efficiency 
in the limit of A —)■ 0. 

The fact that the power vanishes in this limit is the 
same as in the standard heat engine; the Carnot cycle 
produces zero power because the quasistatic limit of the 
Carnot engine takes infinite time for us to operate the 
engine. As A —>■ 0, Jn goes to zero, which means that 
the event that one particle transmits becomes rarer and 
rarer because is the mean value of the particle current. 
In other words, as A —>■ 0, the mean transmission time 
of particles, that is, the average time that it takes for 
one particle to transmit, becomes longer and longer. We 
expect that there is a relation as A ^ ht~^, where t is 
the mean transmission time of particles. It thus takes 
infinite average time for particles to transmit from the 
left reservoir to the right reservoir in the limit A —> 0. 
This physically corresponds to the quasistatic limit of the 
cyclic heat engine, in which it takes infinite time for us 
to operate the engine. 


V. THE CASE OF INELASTIC SCATTERING 

We here extend our theory to the case of inelastic scat¬ 
terings, such as electron-hole and electron-phonon inter¬ 
actions, in which the number of electrons is conserved, 
but their energy is dissipated from the central system. 
In reality, the dissipated energy goes to an environment. 
In order to model this situation theoretically, we intro¬ 
duce a third bosonic reservoir (Fig. i [50] , into which 
the dissipated energy goes. By tuning its temperature Tt 
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and chemical potential /ix, we can make the net particle 
current going into the reservoir zero and the net energy 
current going into the reservoir positive. Note that the 
energy current going into the third reservoir is equal to 
the heat current going into it because the net particle 
current going into it is zero; see the expression of the 
heat current in Eq. Let us denote Qt as the heat 
current going into it hereafter. By additionally tuning 
transmission probabilities going into the third reservoir 
from the other reservoirs, we should be able to set Qt at 
any values observed in experiments. Note, however, that 
this model cannot represent all cases of inelastic scatter¬ 
ings; for example, we may have to consider a non-Markov 
reservoir if the dissipation is non-Markov. We here focus 
on the dissipation which we can represent with the model 
in Fig. We again let Ql and Qr denote the heat flow¬ 
ing per unit time into the left and right reservoirs, re¬ 
spectively, and W the work done by the central system 
per unit time. Note that we define the heat currents as 
Ql = ~Jq and Qr = Jq in the present paper. 

The first law of thermodynamics then gives the follow¬ 
ing relation: 

Ql + Qr + Qt + VT = 0. (25) 

The efficiency is given by 


VI. SECOND LAW, RECIPROCITY, AND 
UNITARITY 

We here remark on the relation of the entropy pro¬ 
duction with the reciprocity and the unitarity. Using 
Onsager’s reciprocal theorem, Benenti et al. [9] recently 
proposed an interesting argument in the linear-response 
regime that a thermoelectric engine could have a finite 
power with the Carnot efficiency under a magnetic field. 
As this proposal is only in the linear-response regime, 
we hereafter argue its possibility in nonlinear-response 
regimes. 

For this purpose, let us first impose Onsager’s recip¬ 
rocal relation without the unitarity of the transmission 
probability. Without the unitarity, however, we can 
choose any expressions of Jn because particles are not 
conserved without the unitarity. For example, suppose 
that we choose seemingly plausible expression of Jn as 

/ OO 

(ie(TL^R(e, B)/L(e)-rR^L(e, B)/R(e)). (30) 

-OO 

Here rL_>.R(e, H) is the transmission probability for the 
electrons from left to right at energy e and TR_,.L(e, i?) 
that for the electrons from right to left at energy e. This 
expression of Jjv leads to the expression of Jq in the form 


W Ql + Qr + Qt (^ , Qr\ , Qt 
—Q h Ql \ Ql j Ql 


■ (26) 


By defining the dissipative heat current as Qt = Jq, we 
can rewrite the efficiency Eq. (261 as 




(27) 


Jq= (ie(e-/XL)(TL^R(e,H)/L(e)-rR^L(e,H)/R(e)). 

J —OO 

... 

Let us further choose the transmission probabilities as 
TL_j.R(e, B) = a{B) at e = Cc and otherwise zero, as well 
as rR_,,L(e, B) — b{B) at e = and otherwise zero, where 
a{B) and b{B) are constants with a{B) > b{B). Note 
that these choices are not prohibited under Onsager’s 
reciprocal relation rL_>.R(e, H) = rR_,.L(e, —H) [47j. With 
Eqs. (301 and (31), we find that the efficiency 


and obtain the work in the form 

W = (-Ql - Qr) - Qt = - J§) - - J^- 

(28) 

In the case of no dissipation Jq = 0, the expressions of 
the efficiency and the work reduce to Eqs. (|^ and ([^ in 
the case of two reservoirs, respectively. Since we consider 
the case in which Qt = > 0, we find that the effect 

of the dissipation decreases the efficiency and the power. 
The entropy production of the whole system is given as 
follows: 


_ ^ Or Qt 
“ Tl Tr + Tt 

7L tR tT 

= _L iQ _L iQ 
Tl Tr Tt ■ 


(29) 


Because Tq > 0, the entropy production is still nonneg¬ 
ative if we introduce the third reservoir. 


IV _ JNjfJ-R — ^J-h) 

jL ~ jL 

•-'Q 


(32) 


becomes the Carnot efficiency rj = rjc, whereas the power 
has a finite value JAr(/rR—/in) = {a{B)—b{B)){p^—jjLi^) > 
0. Therefore, we could obtain a finite power with the 
Carnot efficiency only under Onsager’s reciprocal theo¬ 
rem, which is the nonlinear version of Benenti et al.’s 
argument in the linear-response regime [5]. 

However, the unitarity of the transmission probability 
prohibits this situation. The unitarity condition of the 
transmission probability yields TL_j.R(e, B) = rR_>L(e, B), 
which makes a{B) = b{B) at any B, and hence we find no 
power at the Carnot efficiency. This corresponds to the 
fact that Brandner et al. |14] found a new bound among 
Onsager’s coefficients from the unitarity of the transmis¬ 
sion probability considering the three-terminal model in 
the linear-response regime. This bound prevented the 
power from being finite with the Carnot efficiency. We 
also note that not Onsager’s reciprocal theorem but the 
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unitarity guarantees the non-negativity of the entropy 
production in the mesoscopic transport theory. 

Brandner and Seifert [T3] argued the attainability of 
the Carnot efficiency at a finite power in a multiterminal 
model when the number of probes (a probe is a special 
reservoir in which we set the temperature and chemi¬ 
cal potential so that the net particle and heat currents 
flowing into the probe may be zero) is infinite, but they 
themselves denied it later in Ref. [25] because of a numer¬ 
ically conjectured inequality for the Onsager coefficients. 
From our point of view, we presume that the inequality 
found in Ref. [3S| is based on the unitarity. 

To summarize, Benenti et al. |S] must have missed ad¬ 
ditional constraints that prevent the Carnot efficiency at 
a finite power. In mesoscopic thermoelectric heat engines 
with the Landauer-Biittker formula, the constraint is the 
unitarity of the transmission probability. It is an open 
problem as to what constraints prevent the Carnot effi¬ 
ciency at a finite power in general cases. 


VII. CONCLUSION 


In this paper, we derived the heat currents coming into 
or from the central system under reasonable assumptions, 
with which we can produce Eqs. Q and (14) as those 
of cyclic heat engines. The same expressions have been 
used in the previous researches without detailed expla¬ 
nations. These heat currents gave us interesting results 
when we applied them to the mesoscopic thermoelectric 
heat engine within the Landauer-Biittiker formalism. We 
found that the heat currents in this engine Eq. cor¬ 
rectly give the non-negativity of the entropy production. 
We also found that the only transmission probability to 
achieve the Carnot efficiency is the delta-like function 


Eq. (21) and that the unitarity of the transmission prob¬ 
ability guarantees the Carnot efficiency at zero power. 

It will be interesting to incorporate electron-electron 
interactions to our theory. If the interactions are elastic, 
Eq. (|^ is still valid because it does not break our assump¬ 
tions; we may even be able to derive microscopic expres¬ 
sions extending the Landauer-Biittiker formalism as in 
Refs. [51] [^. Using Christen and Biittiker’s nonlinear 
scattering theory, in which electron-electron interactions 
are treated as mean-field charging effects [TH| |T9| 1351 - 
I37||53H55], we can construct the same formalism by re¬ 
placing T(e) with T(e,rL,'7R,^L,m)- 

It is also interesting to look for stronger bounds than 
the second law of thermodynamics in some systems, 
which make it impossible to reach the Carnot efficiency. 
Eor a system with many probes, for example, especially 
under the broken time-reversal symmetry, the transmis¬ 
sion probability may not have enough degrees of freedom 
to be in the form of Eq. (21). 

We Anally mention the possibility of experimental re¬ 
alization of the mesoscopic heat engine. It may be easy 
to make the setup of the steady-state heat engine, partic¬ 
ularly the mesoscopic thermoelectric heat engine, thanks 
to the improvement of experimental techniques in meso¬ 
scopic transport systems. The results in this paper can 
be used from quantum point contacts to one-dimensional 
nanowires, as well as to cold atoms [55] if we utilize the 
Landauer-Biittiker formula. A possible difficulty is how 
we experimentally observe the heat currents and the ef¬ 
ficiency of the heat engine. This problem may be solved 
in the near future because the technique of observing the 
energy current has been improved recently [57]. 
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